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Abstract 

General two-particle system is considered within the formalism of Fokker-type action integrals. 
It is assumed that the system is invariant with respect to the Aristotle group which is a common 
subgroup of the Galileo and Poincare groups. It is shown that equations of motion of such system 
admit circular orbit solutions. The dynamics of perturbations of these solutions is studied. It 
is described by means of a linear homogeneous set of time-nonlocal equations and is analyzed in 
terms of eigenfrequencies and eigenmodes. The Hamiltonian description of the system is built in 
the almost circular orbit approximation. The Aristotle-invariance of the system is exploit to avoid 
a double count of degrees of freedom and to select physical modes. The quantization procedure 
and a construction of energy spectrum of the system is proposed. 
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I. INTRODUCTION 



Fokker-type action integrals [l, 2| represent an approach to the relativistic particle dynam- 
ics which is alternative or complementary (depending on a point of view) to field-theoretical 
approaches. Known about for a century 3|-[5|, but mainly owing to the Wheeler- Feynman 
electrodynamics , this approach was generalized to other field-type relativistic inter- 



actions including cases of higher-rank tensor fields 8l4lCll|. gravitation llNl4|. confining 
interactions flsUlot etc {20]. 

A variational problem based on the Fokker-type action describes a dynamical system 
with time non- locality, i.e., it leads to difference-differential or integral-differential equations 
of motion for which the Cauchy problem is unsuitable. Consequently, the study the phase 
space (i.e., a set of possible states), the construction of the Hamiltonian description and 
quantization of such a system are non-trivial tasks. 

Serious effort was made to develop Hamiltonization procedure for the Fokker-type action 
integrals. In general, this is attained by means of reformulation of the problem into another 
but time-local form. Here we just mention two such schemes. 

The first is a formal expansion of the Fokker-type action into the Lagrangian action with 



higher derivatives (of order up to infinity) 2l|, |22|, with a subsequent use of a modified 



Hamilton-Ostrogradsky formalism 23|]. In the second scheme developed by Llosa et al. [24, 
251] the variational problem is reformulated into a static one for particle world lines treated 
as temporally extended strings. In practice both schemes can be realized approximately: 
in the first the quasi-relativistic approximations 23| are used, for the second the coupling- 
constant expansion method was developed ^,]. Thus the resulting Hamiltonian description 
of an N-particle system is built on the 6 N- dimensional phase space, as in a non-relativistic 
or free-particle case. 

Among not numerous solutions to Fokker-type variational problems studied in literature 



the class of two-particle circular orbit exact solutions |27H29| is of particular interest. These 
solutions include domain of essentially relativistic motion of strongly coupled particles and 
thus they stand apart the field of application of quasi-relativistic approximations and a 
coupling constant expansion. In the case of the Wheeler- Feynman electrodynamics equations 
for small deviations from circular orbit were derived and studied jsol. The analysis revealed 
bifurcation points in the highly relativistic domain of the phase space where redundant (as 



to compare to a non-relativistic case) unstable degrees of freedo m g et excited. This result 



was approved by direct numerical 



31 



32| and global analytical |33| study of the Wheeler- 
Feynman two-body variational problem. We do not discuss here a physical meaning of 
highly relativistic unstable solutions mentioned above. But on the whole the almost circular 
orbit (AGO) approximation turns out more informative and thus appropriate in the highly 
relativistic domain than the quasi-relativistic or weak coupling approximations. 

Of physical interest is a study, in AGO approximation, of various Fokker-type systems, 
especially those which may have relevance to the relativistic bound state problems in the 
nuclear and hadronic physics. In particular, some Fokker-type systems with confining inter- 
action may serve as relativistic potential model of mesons 15l4l8|. Thereupon the quantiza- 
tion procedure of Fokker-type models which is based on AGO approximation scheme should 
be developed. One can relay in this way on the Bohr quantization of circular orbits {29 1 and 
a heuristic suggestion 3^ to use the Miller's quantum condition [ssl for WKB quantization 
of AGO in the Wheeler-Feynman electrodynamics. 

In this paper we propose a substantiate quantization recipe of a two-body Fokker-type 



problem of general form in AGO approximation. The recipe in based on an implicit E 



amil- 



25|. We 



tonian description of the system which, in turn, is built by means of Llosa scheme 
suppose that a Fokker-type system is invariant with respect to the Aristotle group which is 
a common subgroup of the Galilei- and Poincare groups. By this both non-reativistic as well 
as relativistic systems are involved into consideration. The symmetry with respect to time 
translations and space rotations results in the existence of generalized Noether integrals of 
motions jssl, i.e., the energy and the total angular momentum. We show that, under rather 
general condition, a system admits a circular particle motion with a given constant angular 
velocity. It is taken as zero-order approximation for non-circular motions. In order to study 
a perturbation to circular orbit solution we use a uniformly rotating reference frame where 
circulating particles are motionless. Then AGO solution is formulated in terms of small 
deviations from fixed (and presumably equilibrium) particle positions. We obtain a time- 
nonlocal action principle for these deviations and derive corresponding linear homogeneous 
set of integral-differential equations of motion. Fundamental set of solutions to these equa- 
tions can be expressed in terms of characteristic frequencies and amplitudes of generalized 
normal modes. The amplitudes are shown to be canonical variables, the frequencies are 
functions of the total angular momentum, and all they constitute a correction to zero-order 



circular-orbit energy. Then the quantization is trivial. In general, it must be complemented 
by some selection rules for separation of physical modes out from all variety of them. It is 
discussed in details. 

The paper is organized as follows. In Sec. II we apply the AGO approximation method to 
a single-particle Lagrangian system which is rotary-invariant and local in time. This section 
has a rather methodological meaning since main points of the method are demonstrated, 
and useful definitions and notations are introduced there. In Sec. Ill the method is extended 
for a general Galilei-invariant two-particle system. We consider the latter as a time-local or 
slow-motion limit of a wide class of Fokker-type two-particle systems examined in Sec. IV. 
Perturbations to circular orbits are shown to be described by a linear set of time-nonlocal 
equations of motion. Symmetry and dynamical properties of this set is studied in various 
subsections of the Sec. IV as well as in appendices. In particular, in Appendix D the 
Hamiltonization and quantization of a linear nonlocal system is discussed in detail. 



II. ROTARY-INVARIANT SINGLE-PARTICLE DYNAMICS 

Let us consider a system of single particle which is invariant under the time translations 
t — 7- 1 + A and the space rotations x — > Rx, where A G M, R G 0(3) and x = {xi] i=l, 2, 3} G 
E"^. The Lagrangian function L{x,x), satisfies the equality: 

L{Rx,Rx) = L{x,x) (2.1) 

and thus has the following structure: 

L{x, x) — L{x'^, X ■ X, x^) = L{a, 7) 

Following the Noether theorem the energy E and angular momentum J, 

E^xp-L, (2.3) 
J ^xxp, (2.4) 

are conserved; here p — dL/dx and " X" denotes a vector product. 

The system is exactly integrable with standard methods. We consider this example here 
in order to demonstrate the idea of the approximate method applied below to a Fokker-type 
system. 
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(2.2) 



II. 1. The description in a uniformly rotating reference frame 



First of all we perform the coordinate transformation x ^ z corresponding to transition 
to a uniformly rotating reference frame: 

x{t) = S{t)z{t), with S{t) = exp(tQ) (2.5) 

where Q G o(3) is a constant matrix. We introduce the vector f2 which is dual to Q: 
= ~\^k''^ij- This vector determines the angular velocity Vt = and the direction 
n = Q,/VL oi rotation of a reference frame. 

Using the equality [f^t;]- = (f2 x v)i we complement the coordinate transformation (12. 5p 
by the velocity transformation: 

x = Su = S{z + Q.z) = S{z + nx z) (2.6) 

and calculate the Lagrangian in the rotating reference frame: 

L{z,z-n)^ L{Sz,Su). (2.7) 

This function of z, i; is rotary invariant but with respect to the time-dependent realization 
of 0(3): z — S^^(t)RS(t)z = S(— t)RS(t)z. The corresponding conserved quantity is the 
same vector of angular momentum J as in eq. (12. 4p . 

Besides, the Lagrangian L(z, z\ Vt) does not depend on the time t explicitly, so that the 
conserved quantity E exists although it differs from the energy (12.31) : 

E-^-'i-L. (2.8) 

It is related with the integrals (12. 3p . (12. 4p by means of the equality: 

E = E-n-J. (2.9) 



II. 2. Circular orbit solutions 



Let us consider the solution of the above dynamical problem which is static in the rotating 
reference frame: z = z = R. The Euler-Lagrange equations take the form: 



dl/dz 
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0. 



z=0 



(2.10) 



Taking into account the structure (12. 7p . (12. 2p of the Lagrangian, the equations (I2.10p read: 



RL^ -fix {fix c)L^ = 0, (2.11) 
where La = dL/da e.t.c; in the present case a = R^, (3 = R ■ {ft X R) = 0, and 7 = 

{n X Rf = n^R" - {n ■ Rf. 

Let us consider two cases. In the special case R\\ we have /3 = 7 = while a = B? = 
l-Rp must satisfy the equation aLa{a, 0, 0) = 0. The solution x = SR is truly static: 

X = S(^R^ + g_x_R) = 0. 



In general, j]' f2. Then the eq. (12. lip determines both the direction of the vector R, 

n R = fl x = =^ R±n, x±n, 

as well as the relation of R = \R\ and Q = 

La{R^, 0, n'^R^) + n^L^{R\ 0, n^R^) = O. (2.12) 

Thus R = R{n)R, where R ± ft, \R\ = 1. 

The values of the integrals of motion on the circular orbit solutions are: 

J(°) = 2Rx{nx i2)Lf = 2nR^Lf, (2.13) 
m = (2.14) 
= 2n^R^Lf - = fijW - (2.15) 

they depend on only; here L'^^\ L^^ e.t.c. denote values of corresponding functions on 
the circular orbit solution. 



II. 3. Equations of motion in the linear approximation 

Let us put 

z = R + p, u = v + p + nxp, with V = n X R, 
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where \p\ ^ |-R|, and expand the Lagrangian (12. 7p in the vicinity of the extremal point R 
up to quadratic (with respect to p) terms. One obtains: 

f/ .X w„ .X w„ X dLiR.O) dL(R,0) . 

Liz,z) = L{R + p,p)^L{R,0)+ ' ■p+ \r ' p 

oz oz^ 

Q total derivative 

2 I dz^dz^ dz^dz^ dz^dz^ 

= L(o) + L(2). (2.16) 

(the argument O of L is omitted here). Using the notations 

we write down the second-order Lagrangian 

= i(^uPV'' + 2L,jpV'' + Li^P'p'), (2.18) 
and corresponding equations of motion: 

Li,p^ + (L,j - Lji)p^ - U{fP = 0. (2.19) 

It is convenient to chose unit coordinate orts as follows: €3 ft ^) ^1 Tt -R? ^2 = 
^3 X ei tt f2 X i2, and decompose the vector p = {p^,p'^,p^} into coordinate components. 
Then taking the rotary invariance of the Lagrangian into account (see Appendix A) one 
obtains the equations of motion in the following form: 

Ln/ + (^12 + L^/R)p' - L^ip' - L^ip" = 0, (2.20) 
- (L12 + LilR)p' - Liip' - L^ii? = 0, (2.21) 
-;^(^^V + P^) = 0, (2.22) 

The equation fl2.22p splits out from other ones of this set; it describes the harmonic oscil- 
lations in the direction €3 O with the frequency Vt. Physically, this can be treated (in 
the linear approximation) as a motion of particle along a plane orbit, the normal to which 
differs from €3. Inother words, this mode combines with a circular orbit solution resulting a 
new one with the angular velocity fi = Rfi (where R is a rotation by some small angle). In 
order to avoid double counting of degrees of freedom one can assign the constraint p^ = 0. 



The equation fl2.2ip can be integrated out once: 



(Li2 + L^/R)p' + Lj2p' + L^2P^ = C, 



(2.23) 



with the integration constant C. Let us show that one can put C = without loss of 
generahty. Indeed, if C 7^ 0, the set of equations (12.201) . (12.231) possesses the solution 

= Po, = Pq + Pot with some constants pi, pg which are proportional to C (the constant 
Po falls out the equations and is unimportant). The variable grows beyond all bounds 
of applicability of the linear approximation unless p^ = and p^ = 0. On the other hand, 
the solution with p^ 7^ or/and p^ 7^ can be treated (in the linear approximation) as a 
motion of particle along the circular orbit of the radius R = R + pi with the angular velocity 

= i7 + p1/R, i.e., as some zero-order solution. Thus, it is sufficient to put C = yielding 
P' = Pi 

Apart these two modes (in and €2 directions) which we will refer to as kinematic ones, 
the system (I2.20p - (l2.22p possesses the third dynamical mode which is of physical interest. 
Looking for a bounded solution of the linear set (I2.20p - (l2.22p we use the ansatz 



p' = e\u)e-'^\ 



(2.24) 



real part of which makes a physical sense. Substituting this anzatz into (I2.20p - (l2.22p yields 
the set of equations V{u)e{u) = for a polarization vector e{u) with the dynamical matrix 



where = -^12 + ^2/ 
The secular equation: 
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ILOL 



12 



22 








det V{u) = - Q') {LuL^^ - + [^11^22 - ^id } = 

leads to three solutions for eigenfrequencies squared corresponding to three degrees of free- 
dom of the system. The only one solution corresponds to the dynamical mode: 



2 _ ^12 ^ ^11^ 



22 



-^ii-^22 ~ -^i2 



(2.25) 



If uf > 0, both values of eigenfrequencies ±ui are real; they permit us, using the ansatz 
fl2.24p . to construct a physically meaningful real and bounded solution of equations fl2.20p - 

Two other eigenfrequencies squared and corresponding eigenvectors, 

Loj = 0, e = {0,1,0}; 

ul = Q\ e = {0,0,1}, 

are images of the kinematical modes mentioned above. 



II. 4. Integrals of motion in the linear approximation 

Let us start from the angular momentum J. It can be presented as follows: 

f)f 

J = SzxS— = Sr, (2.26) 
oz 

where the vector 

is not, in general, conserved. Its components, in the linear approximation, can be presented 
as follows: 

= + p^'KL^ + L^,p^ + L^J + ...) 

+ tI^\ (2.28) 



The components Tf^'' and T^^'' have the explicit form: 

j{0) ^ j{0) ^ ^ ^^^^ ^2.29) 

= L^p\ = -^L^p', (2.30) 

t!^'^ = ^ {L^2P' + Lap' + ^22P'} = ^- (2.31) 

It is evidently that J*^'^^ = ST'-^-* = T'-^-*. Besides, the only kinematic modes but not 
the dynamical one contribute in T*-^-* and thus in J*-^^ = S T'-^^ It was pointed out in the 
previous subsection that we can put = and C = without loss of generality. Then 
j-(i) _ 5 2"(^) = 0, and, in the given approximation, J ^ J'^^\Q,), where the function J^^\Q,) 
is defined implicitly by fl212|) -f l213|l . 



Now we consider the energy of the system. First of all, we calculate the correction to the 
zero- order term fl214D of the integral E (12:81) : 



E ^ + ^(2) = + H^yPV' - L,,p'p^}. (2.32) 

So that the first nontrivial correction E^'^'^ to E^^^ is quadratic in p*. It is evidently conserved 
by virtue of the equations of motion in the first-order approximation fl2.20p -( l2.22p . 
Further we are interested not in the integral E but in the energy 

E = n-j + E = z-^ + n-^-L^n-j- + e^^\ (2.33) 

OZ Oil 

It follows from the equalities f l2:28|) -( l232|) : 

E^'^ = IWp' + L,f^P'} + (2-34) 

2 J3 

where i, j = 1, 2. On the other hand, within the given accuracy 

I 2J J 24"^ 



Thus we obtain a useful equality 

E^1]J-LW+E(2), where E^'^ = E^'^ 
which holds with accuracy up to quadratic terms in p's. 



33=0 



(2.35) 



II. 5. Hamiltonian description and quantization. 

The Legendre transformation p 1— )■ tt = dL/dz leads to the Hamiltonian description 
with the Hamiltonian function H{p, tt; f2) to be the integral of motion E (12. 8 p in terms of 
canonical variables p, tt. 

The fixed auxiliary vector O was introduced to specify the rotating reference frame and 
then the circular orbit solution. In order to generate a set of all possible circular orbit 
solutions we let f2 to be a variable of angular velocity. It follows from fl2.27p and fl2.33p that 
the Legendre transformation with respect to both z = p and f2 leads to the Hamiltonian 
description with the Hamiltonian function H{p, tt; Y) to be a conventional energy. Rotary 
invariance and a Hamiltonian constraint born from the identity (I2.27P provides a proper 
balance of degrees of freedom in the phase space enlarged with the T variable. 
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It is convenient in our case to proceed from the expression fl2.33p . First two terms depend 
on J = I J| = \ T\ and VL only. Using f l2.12p . fl2.13p one can express Vt = f{J). Thus, in 
zero-order approximation, we have the Hamiltonian: 

H^'\j) = E^'\n,J)\^^^^j^. (2.36) 

Similarly, coefficients Ly, of the quadratic form E^"^^ in fl2.35p turns into functions of 
\J\. We note that within the Hamiltonian description the components Jj of the angular 
momentum J satisfy the Poisson bracket relations (PER): 

{Jii Jj} ^ij Jk 1 

non-triviality of these PER is due to the fact that original variables are not velocities but 
quasi- velocities. 

To complete the Hamiltonization one should eliminate in i?'-^^ velocities (i = 1, 2) in 
favour of canonical momenta tti = dU-'^^ jdp^ satisfying the PER {pSttj} = 5] (others are 
trivial). For a quantization purpose it is better to use normal mode complex amplitudes 
satisfying PER: {AcA^} = — i^Q,^ (a = 1,2 in our case). Then, using results of Subsec. 
II. 3-4, the energy correction is put in the Hamiltonian form: 

1/(2) = UJ^\Aj^ = UJr{J)\Ar\'^. (2.37) 

' ^ a 

Here ujr = cui (>0) is the characteristic frequency of dynamical mode (12.250 : we redenoted 
the subscript a = 1 — t- r to hint that this degree of freedom corresponds to radial oscillations. 
The radial frequency uJr{J) is expressed in terms of J rather than Q. 

We note that the kinematic mode with the frequency a;2 = is not oscillator-like and 
must be suppressed; the corresponding solution is discussed in Subsec. II. 3 after eq. (12.230 . 
The corresponding contribution in the Hamiltonian (I2.37P drops out automatically. 

Finally we have the Hamiltonian H = H^^^ + H^"^^ which is ready for quantization: 
variables are replaced by operators and then - by their eigenvalues as follows: 

J — )• J; Ar ^ 

J ^ ^ V^£(7Tl) , £ = 0,1,...; 

\Ar\'^ ^ \{Aral + alAr) ^ rir + \, n^ = 0, 1,... (2.38) 

It is implied, due to a perturbation procedure, the condition H^"^^ <^ H^'^^ which is mainly 
satisfied by < £. Then ^/e(JTT) ^ £ + ^. 
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General structure of the Hamiltonian fl2.36p . f l2.37p and its spectrum agree completely 
with corresponding results derived with the standard Hamilton- Jacobi and WKB methods. 
In Appendix B the quantization method is demonstrated on the example of a nonrelativistic 
particle in a power-law potential. 

III. GALILEI-INVARIANT TWO-PARTICLE DYNAMICS 

The rotary-invariant single-particle dynamics studied in the previous section gives us an 
important tool for the description of a two-particle system. Any isolated system, anyway 
the non-relativist ic or relativistic one, possesses 10 conserved quantities: the energy E, the 
momentum P, the angular momentum J and the boost K. This is consequence of invariance 
under the action of a symmetry group: the Galilei group in the non-relativistic case, and 
the Poincare group in the relativistic case. It is known, both in the non-relativistic and 
relativistic cases, that in the rest reference frame fixed by the condition P = (and also 
1^ = in the classical, i.e., non-quantum description) the two-particle dynamics can be 
reduced to an effective single-particle one with a residual symmetry group to be 0(3) xT, 
(here T denotes the time translation group). Below we consider a non-relativistic two- 
particle system and reduce it to an effective single-particle one in the AGO approximation. 

III.l. General dynamics and circular orbits 

The Galilei-invariant two-particle Lagrangian has the following general form: 

2 2 

L{xi, X2, xi, ±2) = -Y^l + F{x^, X ■ X, x^) = Y ~Y^"^ ^ ^^'^^ 

a=l a=l 

where x = Xi — X2- The corresponding 10 integrals of motion are: 

Uln 9 OF 



2 " dx 

a=X 
2 



P = YmaXa, (3.3) 



2 



J = ^ niaXa XXa + XX. — , (3.4) 



dx 

2 



K = YmaXa-tP. (3.5) 



a=l 
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= 0, a = 1,2 

i,=0 



Non-inertial variables are introduced similarly to the single-particle case: 

Xa = SZa, Xa = SUa = S{Za + 0^ >^ Za) =S{Za + Va). (3.6) 

In these terms the Lagrangian 

Z2, zi, Z2] O) = L{z, ui, U2) (3.7) 

does not depend on time t explicitly and thus it generates the corresponding integral of 
motion: 

^ = E^'^-^-^' (3-8) 

related to the original integrals ( 13. 2p . (13. 4 p by means of eq. (12. 9p . 
Circular orbit solutions are determined by the conditions: 

dl/dz^ 

which explicit form is: 

-m^fi X i;a + 2(-)^(2F„-ri X vF^) =0, a = 3 - a. (3.9) 

Multiplying left- and right-hand sides of these equations by f2 yields: 

2(-)^0 X zF^ = {) =^ z ±n. (3.10) 

We note that eqs. (13. 9p are invariant under translations along O, i.e., under the trans- 
formations z'^ = Za + An with an arbitrary A G M. Indeed, it is evidently that z z' = z, 
and also 

V'^ = n X Z'^ = n X {Za + \n) = n X Za = Va. 

Taking into account the equality (I3.10p one finds: 

n- Za = ^-{Za-{Za- Za)) = O ■ Z^, 

i.e., zl = n ■ Za = zl but no information for the last quantity follows from eqs. (13. 9p . Thus 
one can choose 

n-Za = o =^ za±n, 

which simplifies the system (13. 9 p to the form: 

{rrnn^ + 2[F^ + n^F^]}zi - 2[Fa + n^F^]z2 = 0, 
-2[F„ + n^F^]zi + {msfi' + 2[F„ + n^F^]} Z2 = 0, (3.11) 
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where a = z^, (3 = and 7 = Q'^z'^. 

The hnear homogenous set of equations (13. lip possesses a non-trivial solution if its de- 
terminant vanishes: 



n'^ {mimsl^^ + 2(mi + m2)[F„ + fi^F^]} = 0. 



(3.12) 



This is a relation between Q and z. One solution of eq. (I3.12p is: i7 = 0. In this case the 
set (13.91) reduces to the equation F^z = 0. If an interaction is not singular at 2 = then we 
have the solution 2 = 0. Otherwise, the solution is determined by the equality = which 
is the condition of extremum (here, the minimum) of the static potential of interaction. As 
in the single-particle case, this solution is not suitable. 

Other roots of secular equation are determined by the condition: 



(here /i = 17111712/ {mi + 1712) is the reduced mass) which being combining with (13. lip yields 
the set: 



The relation of R = \R\ and Q is, evidently, defined by (I3.13p . 

The values of the integral of motions (I3.2p -( l33|) on a the circular-orbit solutions are as 
follows 



they obviously correspond to a rest of the system as a whole. 

III. 2. The dynamics in the AGO approximation 

Similarly to the single-particle case one puts: 

Za = {-fRa + Pa, Ua = Va + pa + ^ X Pa^ wherC Va = (-)"0 X Ra, 

and expands the Lagrangian (13.71) in p^, p^. One gets L ^ L^^^ + L^^) where 



fj,n^ + 2[F„ + n^F^] = 0, 



(3.13) 




(3.14) 



J(°) = nn^if, + 2Ff ), = R^n^lf, + 2Ff ) - F^; 



(3.15) 



ab 
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„2f (0) 

with the coefficients L„,- hi = tt-tt 



-laibj = r^^i etc. The corresponding equations of motion, 

^ ] {Laibj pI + {Laibj '~ Lbjai)Pb "~ LaibjPb) ~ 0' 
b 

have the following explicit form: 



where F, = j^^^ 



(0) 



etc., = Fjj — Fji and p = — pg- Summing up these equations 
over a = 1, 2, ffist with the wight 1, and then with (— )°' ^^"^°^_^ , splits these equations: 

- Qi^jQ^ + 2e,j''VtkQ' -Qi = Q (3.16) 

[^^{^f5,, - n,nj) + Fij] + + %]]p^' - [fi5,j + FijXfP = o, (3.17) 

where o = V„ — — p„ is a deviation of the center-of-mass position. 

Let us consider the equation f l3.16p . Choosing orts as in the 1-particle case simplifies it 
to the set: 

n^g^ + 2VLq^ - = 
-g^ = 0. 

In order to cut unbounded solutions off we search a solution in the form: g'^ = 5*e~"^*, and 
arrive at the secular equation: u'^{u)'^ — Q'^)'^ = 0. We claim without details that eigenvectors 
belonging to the degenerate eigenvalue u"^ = correspond to a rotation of the vector g 
with the frequency Q. In the fixed (motionless) reference frame this solution is the constant 
vector S-Lfl. The u = mode possesses constant eigenvector e\\fl. All three modes can 
be compensated by the translation of the origin of coordinates, i.e., by redefinition of the 
center-of-mass reference frame. Thus these modes are kinematic, and one can put ^ = 0. 

The set of equations fl3.17p can be obtained from the set (12.191) or (I2.20I) - (I2.22I) by means 
of formal substitution L — >■ L, where 

L(a,A7) = i/i7 + F(a,/3,7). (3.18) 

This set leads to one dynamical mode with frequency fl2.25p (with the change L L) and 
two kinematic modes in addition to three ones described just above. Particle eigenvectors 
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Cq of all the kinematic modes have the following components: 



W2 = 0, ei = {0,i?i,0}, 62 = {0,-i?2,0}; (3.19) 

co3 = ±n, ei = {0,0,i?i}, 62 = {0,0,-i?2}; (3.20) 

co4,5 = ±^, 6i = {l,Ti,0}, e2 = {l,Ti,0}; (3.21) 

006 = 0, 61 = {0,0,1}, 62 = {0,0,1}, (3.22) 

where Ra = \Ra\ and Ra (a = 1,2) are defined in (13.141) . 

After the kinematic modes are suppressed, the subsequent analysis is reduced to the 
single-particle case considered in Sec. II with the effective centre-of-mass Lagrangian fl3.18p . 
as it is in the standard treatment. 



IV. TWO-PARTICLE FOKKER-TYPE DYNAMICS 



In this section we extend the AGO approximation method to the formalism of action 
integrals of Fokker type. We start with a two-particle Fokker-type action of general form 
38|: 

2 



dtaLa{ta,Xa{ta),Xa{ta)) 

a=l ^ 

+ jjdtidt2^ {ti, t2, Xi{ti), X2{t2), ii(ti), X2{t2)) . 

The variational problem leads to the integral-differential equations of motion: 

d d d 



(4.1) 



dZa dta dZa 



(La + A,) = 0, a = 1,2, 



(4.2) 



, d d . d .. d , 
where -— = — h Xa ■ t: h Xa ■ and 



dt„ 



dta 



dXa 

Ai 



dXa 

dt2$, 



Ao 



dti $. 



(4.3) 



For a physical reason we are interested mainly in the case where the system is invari- 
ant under transformations of the Aristotle group [36| (see also (37|), i.e., time and space 
translations and inversions as well as space rotations. The Aristotle group is a common 
subgroup of the Galilei and Poincare groups. Thus this case includes both non-relativistic 
and relativistic non-local systems into consideration. 
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IV. 1. Symmetries and conserved quantities 



Symmetry properties of the action (14.11) determines general structure of the functions 
La(ta,Xa,Xa) (farther referred to as the 1-Fokkerians) and ^{ti,t2, Xi, X2, Xi, ±2) (referred 
to as the 2-Fokkerian) and leads to an existence of integrals of motion studied, for non-local 



[i.e., Fokker-type) systems, in 



The invariance under time translations, t — )■ t + Aq (Aq G M), results in the conditions: 



A-„^i.-g^ = 0, (4.4) 

^«'^ = ilw^° ^ *(«i,«2,...) = *(«i-«2,...) = 4(!',-) (4.5) 

a=l 

and yields the energy integral of motion: 

^ = E ■ - 1 } + ^'^) +ff^t^ (4.6) 



a=l 

#i'ti pca POO pt2 

^ I ■ (^-^^ 

J —00 J t2 J t\ J —CO 

The invariance under space translations^ x ^ x + \ [X ^ M^), yields the conditions: 
X'^U = ^ = 0; (4.8) 

OXa 

^^^ = Y.^ = ^ =^ H-,Xl,X2,-)=H-,X,-X2,..) = H..,X,..) (4.9) 

a=l 

and the conserved total momentum: 

^ = E ^ + - H^t^ dt, (4.10) 

The rotary invariance, 



, OXa Jf dx 

a=l 



La{RXa) = LaiXa), (4.11) 

$(^?, Rx, Rxi, RX2) = $(?9, X, xi, ±2), (4.12) 



yields the infinitesimal conditions: 



X^La^±aX^ = 0, (4.13) 



a=l 
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and results in a conservation of the angular momentum of the system: 

a=l 

- iffdt, d(,{(x, + a.,) X A + X A _ ,,x^} (4.15) 

The consequences of discrete symmetries with respect to space inversions and a time reversal 
will be considered farther. 



IV. 2. Fokker-type dynamics in a uniformly rotating reference frame 

Using the non-inertial change of variables: 

Xa{ta) = S{ta)Za(ta) = SaZa{ta), Xa = Sa^a, (4.16) 

where S{t) and u are defined in (12.51) . (13.61) . and symmetry properties (I4.4p - (l4.5p . (I4.8p - (I4.9I) . 
dlTU-dlH]) one can define "tilded" Fokkerians: 

La{Xa) = LaiSaUa) = La{Ua) = La{Za, Z^, ^) , (4.17) 
$(^9, X1-X2, ±1, ±2) = $(^9, SiZi-S22;2, Silti, S2W2) 
= $(??, SJS1Z1-Z2, SjSini, U2) 
= $(^9,S(^)2l-Z2,S(^)lil,n2) 

= ^^,Zi,Z2,Zi,i2;n). (4.18) 

It is obviously that "tilded" Fokkerians are invariant under time translation. Thus the 
corresponding integral of motion exists: 

^ = E {i- I: - 1} (£« + A„) H-|di, di, A* (4.19) 

where the relations of and $ are similar to (14. 3p . The equality (12.91) holds in the present 
case too which fact can be examined directly with the use of eqs. (14.61) . (I4.15P and (I4.16p . 



IV. 3. Circular orbit solutions 

Proposition. If Fokkerians are invariant with respect to the action of the Aristotle group, 
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i.e., the equahties f Oj) . fH3|) . f Oj) . flOj) . fHTTSj) . f HTTil) hold, the corresponding equations 
of motion f l4.2p posses a circular orbit solution with characteristics described below. 

Proof. Fokker-type equations of circular orbit motion (i.e., equations of a rest in terms of 
variables z^) have the form: 

d 



dZa 



La + K 



0, 



a = 1,2, 



(4.20) 



where L„ = Z,(z„, 0; n) = L(°)(z,; fJ), a = 1,2, 

~K = kt-a $(ti-t2, ^1, ^2, 0, 0; rj) (a = 3 - a) 

= zi, Z2, 0, 0; O) = AW(zi, Z2; O), 

and the last function is common for both values of a = 1, 2. Thus, equations of a rest f l4.20p 
can be presented in the effective Lagrangian form: 



^ ^ I y + A(°) 1 =0, 



a = 1,2. 



(4.21) 



va=l 



In order to take into account symmetric properties of Fokkerians it is convenient to 
represent their general functional form as follows: 

La{Xa) = La{xl) = La{-fa), a = 1,2, (4.22) 

$(■(?, X, xi, X2) = ^{'d, x"^, X ■ xi,x ■ X2, x\, x\, Xi ■ ±2) 

= <l>(t9,a,/3i,/32,7i,72,5). (4.23) 
Invariance with respect to time reversal causes the property 

$(-^,...,-/3i,-/32,...) = *(^,-,/3i,/32,...); 

the space parity is provided automatically. 

In Appendix C scalar arguments a . . .6 are expressed in terms of non-inertial variables 
Za and their derivatives. Using flCTTll -flClil) in (141221) -f H:23D then yields "tilded" Fokkerians. 

For a circular orbit problem it is sufficient to consider the static case z^ = 0. Then the 
expression fIC.ll) for a does not change while remaining scalars flC.2l) - flC.4p simplify: 

= /jf) = Q[zi. sin(Oi?) - (n, Zi, Z2) cos{Q^)] , 



Ta 



5(0) = n'zi 



■ Zn 
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Upon these equalities the 2-Fokkerian takes the following general structure: 



$(0) ^ $(0)^^^ ^2^ |^±|2^ |^±|2^ ^± . ^± (4^24) 

Upon accounting the temporal reversability, 

= (4.25) 
the integrating this function over d yields a general structure of A'^^); 

Since (n, Zi, Z2) = ±|2i X22I, so (n, Zi, ^2)^ = |^i"Pl^2'P ~ i^i ' ^2)'^-. ^^e final structures 
of the 2-Fokkerian and then of the effective Lagrangian involved in eq. fl4.2ip is expressed 
via four scalar arguments: 

A(°) = AW(z2, \z{\\ |4|^z^4;^])^A(°)Kal,a2,a3;^), (4.26) 

2 

L(°) = L(°)(ao, ai, a^, as; = L^!\o,- Q) + AW(ao . . . ^3; Q). (4.27) 

a=l 

This form is useful for a study of circular-orbit equations (14.211) . 

Using the notations ki = dL^^^ /dai (z = . . . 3) brings (I4.2ip to the form: 

2{-ykoZ + 2kazi + k3Z^ = 0, a = 1,2. (4.28) 

Scalar product of these equations with n yields the condition: 

kon ■ z = =^ n ■ z = =^ z±n 

(we do not consider the solution ko = 0; roughly it corresponds to a static case in a conven- 
tional meaning). 

Notice that z and thus ctq is translation-invariant; zf and thus o"a, cts are invariant under 
translations along n. Thus the solution of eq. (14.281) is specified up to arbitrary vector along 
n. We fix it by means of the conditions: 

n-Za = 0, a = 1,2 =^ Za-Ln, Za = z^. 

From now on both vectors Za are placed on a common plane and the superscript "_L" can 
be omitted. 
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The rest equations ( I4.28P take the form: 

2{ko + h)z, - {2ko - h)z2 = 0, (4.29) 

-{2ko-h)zi + 2{ko + k2)z2 = 0; (4.30) 

they possess non-trivial solution provided: 

4(fco + ki)iko + k2) - (2A;o - k^f = 0. (4.31) 

Then Zi\\z2. Choosing orts of a moving reference frame as follows: €3 = n H 17, ei || Zi, 
€2 = 63X61, one can recast (I4.29l) - (l4.30p into equalities: 

Zi = RiSi, zi = -R2ei, (4.32) 
R2 ^ ks - 2ko ^ 2{ko + k2) 

Ri 2{ko + ki) ^3 - 2fco ■ ^ ' ' 

By f l4.32p we presuppose Z2 ti ^1 ^'Hd thus -Ri > 0, -R2 > 0, as in the Galilei-invariant 
two-particle case. Otherwise i?2 < but this case is rather nonphysical. 

Relations f l4.3ip and f l4.33p form the set of equations determining Ri and R2 as functions 
of ^] ■ 



IV. 4. Integrals of motion along circular orbits 



Static character of circular orbit solutions implies that Fokkerians L^a^ and A^^^ depend 
on the constant vectors (14.320 . Besides, ^^^^ is a function of 'd = ti — t2. Thus integrals 
of motion can be evaluated explicitly. At that it is useful the following "skew" integration 
(14. 7p ansatz valid for arbitrary function /("i?): 

r>00 I't2 



dh dt2f{^) 



t-] POO 



00 Jt2 



d^' {^-^')f{^'). 



Taking into account the time reversal (I4.25P of ^^^^ yields easily an evaluation of the 
integral E^°^: 

2 fOO o 

^(0) = -J2 (4°) + AW) - / d^ 

a=l 

= _ ^ - 2A(°) + d^ = ~Y1 - ^^'^ = -^^'^ ■ 

a=l a=l 



(4.34) 
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An evaluation of the angular momentum is cumbersome: J^^-* = nJ*^°^ where 
J^^'> = 2nY,Rl{La + K)f^ - R1R2 I d^{[(2<l>,-fi2$,)^ + $^^+$^^]sin(fi^9) 

+ [2$5 + {^p, + <l>^J?9]ficos(fi^9)}(°); (4.35) 

the subscripts a, denote derivatives, $0, = d^/da etc. 
It follows from (I4.34p and f l4.35p the following relation: 

Besides, the relation (12. 9p gives the possibility to evaluate the energy: 



^(0) = ^(0) + ^j(O) ^ _^(0) ^ ^5^{0)/5^, 

The linear momentum integral vanishes: P^^^ = 0. 



(4.36) 



IV. 5. Equations of motion in oscillator approximation 



Small perturbations Paita) to circular orbits are introduced naturally: 



Za{ta) = {-TRa + PaiQ, ^a(ta) = PaiQ, a =1,2 



and then substituted into the action (14.11) . Expanding the Fokkerians up to quadratic terms 

in paita), paita) ylclds: 



I = Y,jdtaL^^^ +jjdt,dt2^ 



V 



(0) 



■Pa + 



dZn . " dZa 



(0) 



■Pa 



total derivative 



/(2) 



+ \Y1 J"^^^ {LaijP\p>a + 2Laijpipi + LaijP\p>a) 
a 



with the coefficients Lain = tt-tt 



(0) 



dzidzl 



(0) 



etc. 



(4.37) 
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The equations of motion have the form: 



£„,,pi(t) + - CaiifPAt) + j^t' - t')pi{t') = 0, (4.38) 



where Ca[ij] = Caij - Caji {a = 1,2, a = 3 - a), 



mj Ldij -\- A.aij ) ^aij • • • ; 



(4.39) 
(4.40) 



and $('(9) = d^/d'd. Due to the time reversabihty of the dynamics, the kernel possesses the 
properties: 

Putting Pa{t) = e^e~"^* leads to the set of equations: 



with the 6x6 dynamical matrix 
Dioo) -- 



'^Daibj{uj)ei = 

6=1 



(4.41) 



where the off-diagonal entries Eaji{uj) = f d-i? Saji('i9)e"^'' possess the properties: 



(4.42) 



.CO 



at] 



(4.43) 



The equation (14.411) determines characteristic modes of the system. In particular, eigen- 
frequencies are derived from the secular equation 

det D{uj) = 0. 

Subsequent description of the system depends considerably on properties of the dynamical 
matrix (14.421) . 

First of all we note that the equations of motion (14.381) are not ordinary 2nd-order dif- 
ferential set but they form an integral-differential set which complicates to a great extent 
the analysis of the dynamics. In particular, the Cauchy problem is not appropriate and the 
Hamiltonization is not straightforward. On the other hand, the set is linear which, in turns, 
simplifies somewhat the analysis. The Hamiltonization scheme of a general linear system 
defined by a non-local action integral is discussed in Appendix D. 
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IV. 6. Symmetry properties of the dynamical matrix 



Similarly to the cases of ordinary single- and two-particle systems we should separate the 
dynamical and kinematic modes of the dynamical matrix. This can be done by taking the 
Aristotle-invariance of the system into account. 



Proposition. The dynamical matrix (14.421) built with the Aristotle-invariant Fokkerians ad- 
mits the eigenfrequencies and eigenvectors (I3.19p - (l3.22p . where Ra {a = 1,2) are determined 
by the equations flOT]) . fOS]) . 



Proof. Multiplying l.-h.s. of equalities (gS]), dUD, (1413|) and ( KW) by S^^ = Sj = S'^(tJ 
and expressing them in terms of noninertial variables yields the equalities: 



^ Ink ^ I r 



0, 



vRf _ ^ I j 



d 



d 



d 



~k ^ I Qk n 1 



d_ .k_d_ 

dzi ^"""dzi 



d 



dzl 



a 



VP ' </ r*^ U 7k L Clkp n rr 



dzi 



(4.44) 
(4.45) 



(4.46) 



(4.47) 



here arguments Za-, Za of La and $ are omitted for brevity. 

Symmetry conditions fl4.44p - fl4.45p for 1-Fokkerians result in useful consequences: 



(0) 



(0) 



0, 



0, 



d_ 

dzi 
d_ 
dzi 



1(0) 



(0) 



0, 



0, 



d_ 

dzi 
d_ 
dzi 



liO) 



(0) 



0, 



0, 



where the superscript "(0)" means "on circular orbit solution", i.e., taking the conditions 
zi = Ri, Z2 = —R2, Za = into account. These equalities impose constraints f]E.ip - f]E.6P 
for the quantities Lai, Laij etc. shown in Appendix E. 
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For the 2-Fokkerian we are interested in the following consequences of fl4.46p . fl4.47p : 



i(o) 



0, / 



ozi 



(0) 



-oo 

oo 



d^ 



(0) 







d^ 



d 



dzi 



(0) 



which impose constraints for the quantities Kai = $aj(0), Aaij = $ajj(0), . . .and $aiaj(±0) 
etc. Using the explicit form for the matrix S{^): 

Sj{^) = cos{m)Si + {1 - cos{m)}nin^ - sm{m)n''ej , 

and taking into account the equality S^^(7?) = S"'^('i9) = S(— and the fact that the Fourier- 
transform ^{u) of $(-(9) and its derivatives ^ai{^) etc. are pair-vise functions of u, one can 
arrive at the equalities flE.7l) - flE.12l) in Appendix E. They lead together with the equations 
(]E.ip - (]E.6|) and the equations of a rest Cai = to the following linear relations for elements 
of the dynamical matrix D(0) and D(f2): 



(4.48) 
(4.49) 
(4.50) 
(4.51) 



(4.52) 



^ai3 ~l~ Aaia3 0; 
Ra^ai2 ~ Rd-^aia2 = 0, 

i?,£,[,3] - R,%,,^]m = 0, a = 1, 2, t = 1,2,3 

- ^^3^'%^-ak]i^) + ^'^a^^m =0, j = 1, 2. 

The relations ( I4.48P provide the existence of the eigenfrequency and the eigenvector (13.220 . 
the relations fICTD - of ( 13191) . the relations (H30l) -( H3T1) - of ( Km . and the relations (1432]) 
- of (I32ID ■ 

Thanks to this proposition, the kinematic modes separate from physical modes, eigenfre- 
quencies of which can be found from the equation: 

det D(u) 



Since D(co') is a 6x6 matrix, entries of which are not, in general, polynomial, the secular 
equation (14.530 is rather comphcated. It can be simplified somewhat due to the following 
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Proposition. The following identities hold: 

-Dajfe3(w) = Dasbjiuj) = for any a, 6 = 1, 2 and j = 1, 2. 



Proof can be completed directly using the representation fl4.22p . (14.231] of 1- and 2- 
Fokkerians and examining the following equalities: 



da 



(0) 



d5 



dzl 



d'^a 



dzidzl 



(0) 



(0) 



0, 



0, 



da 



dz^ 



(0) 



d5 
dzl 



(0) 



0, 



dH 



(0) 



0, 



j = l,2, 



dzidzl 

where a, . . .5 are defined by ( lC.l[) -( ICl4l) in Appendix C, and the superscript "(0)" denotes 
the value of a marked quantity on the circular orbit solution ■ 

Thanks to this proposition, det D(co') splits into two factors: 

det D(u;) = det D^(w) ■ det Dll(a;), 

where D^{uj) = \\Daib;{u)\\ (i,j = l,2), = \\DaMuj)\\ . 

The 2x2 submatrix Dll(a;) possesses two kinematic modes (13.201) and (13.221) while the 4x4 
submatrix D-^{uj) - another three kinematic modes (13.191) . (I3.2ip and one the dynamical 
mode. The frequency of the latter can be determined frOm the reduced secular equation: 

det D^iuj) 



u;2(a;2_Q2)2 =0- 

The secular equations (I4.53P or (14.541) simphfy in the case of equal particles. 



(4.54) 



IV. 7. The dynamics of the equal particle system 

The equal particles system is defined by Fokkerians of the following properties: 

La{Xa) = L{Xa), (4.55) 

$(t9, Xi, X2, Xi, X2) = $(-??, X2, Xi, X2, Xi). (4.56) 



Proposition. Along with a plausible assumption the equations of motion (I4.20p for equal 
particles possess a circular motion solution of the form: 



Zi = R, Z2 = —R (and Zi = Z2 = 0) 



(4.57) 



26 



with characteristics described below. 



Proof. It follows from f l4.57p and f l4.29p . ( ]4.30p the equality: ki = k2 which turns into 
identity provided fl4.57p holds. Then r.h.s of the equation f l4.3ip for \R\ = R{Q) factorizes: 

{Ako + 2ki -k^) = Q or/and {2ki + ^3) = 0. (4.58) 

Although solutions of both the equations fl4.58p must be examined, the only first equation 
seems to make a physical sense since it includes the derivative fco = dL^^^ /dz^ which is 
intuitively related to a fource of interparticle interaction ■ 

Proposition. If the equal particle system is invariant under the space inversions the entries 
of the dynamical matrix satisfy the equalities: 

D2i 2j (w) = Dii ij (w) , D2i ij (w) = Dii 2j (u) (4.59) 
Proof. If 1-Fokkerian is invariant under space inversions, 

L{-z,-z) = L{z,z), 
then its second-order derivatives are invariant too. Thus we have 

= Mji-R^ 0) = 0) = (4.60) 

and, similarly, 

L2ij = Liij, L2ij = Liij. (4.61) 
The 2-Fokkerian for equal particles can be presented in the form 

$(??, Zi, Z2, Zi, Z2) = ^F(i9, Zi, Z2, Zl, Z2) + ^F{-19, Z2, Zi, Z2, zi); 

here F{i),x,y,u,v) is some function (for example, the ^{^,x,y,u,v) itself) which is 
inversion- invariant , 

F(i9, -X, -y, -u, -v) = F(?9, x, y, u, v), 
and thus its second-order derivatives are so. We have the equality 

<f 2.2,(^9) = 2w(^, R, -R, 0, 0) + 0' 0) 

= Shi^^ -R^ 0, 0) + 2mi-^, R, R, 0, 0) = $i.i,(-^) 
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and similar equalities for other derivatives and their Fourier-transforms: 



(4.62) 



Then using the definitions fO^ . fICTD . (102]) and the properties fl05D . flOTD . fOT]) and 
fl462|) resuhs in the equalities fH759|) ■ 

Thanks to equalities fl4.59p the set of equations fl4.4ip splits into to subsets: 



It is easy to verify that the subset (14.631) possesses three kinematic eigenfrequencies and 
eigenvectors (I3.2ip . (I3.22p while the subset (I4.64p - another two kinematic modes (13.190 . 
(I3.20p and one the dynamical mode. The frequency of the latter can be determined from 
the reduced secular equation: 



is the reduced dynamical matrix. In the next subsection we discuss possible solutions of this 
equation as well as of more general equations (I4.53P or (I4.54p . 

IV. 8. Predictive treatment of the Fokker-type system. 

Let us call a particle system as predictive if it possesses three degrees of freedom per 
particle. For example it is the Galilei-invariant Lagrangian two-particle system considered in 
Sec. III. Consequently, in AGO approximation the corresponding dynamical matrix admits 
6 modes with real frequencies, 5 of which have a kinematic origin (they are or ±fi), 
and only one mode characterizes a specific dynamics (it corresponds to radial interparticle 
oscillations) . 

The same kinematic modes arise in Aristotle-invariant Fokker-type 2-particle system and 
thus, in a Lagrangian system (as a particular case when 2-Fokkerian includes as well 

as in the Poincare-invariant Fokker-type system (as a particular case with extra Lorentz- 
invariance). 



Dij{u)e^ = with = Cuj - iu;£i[ij] + w^Ajj + ^^^(w) 

= Diiij{uj) + Dii2j{uj), 6^ = ^(e{+e^); 

Vij{u)e^ = 0, with Vij{uj) = Cuj - i + u^Cuj - Hhj(w) 

= Diiij{uj) - Dii2jiuj), = e\ - e\. 



(4.63) 



(4.64) 



det V^{uS) 



where V^{uS) = ||Aj(w)|| (i,j = 1,2) 
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In contrast to predictive systems, the Fokker-type dynamical system may possess an infi- 
nite number of degrees of freedom, due to a time nonlocality of equations of motion. In AGO 
approximation the quadratic term I^"^^ of the action fl4.37p (i.e., 3rd and 4th strings) is of 
the Fokker-type too. Consequently, the dynamical matrix of such systems is not polynomial 
and may possess an infinitely large number of modes with real or/and complex frequencies. 
It has been shown in Subsec. IV. 6 and IV. 7 how kinematic modes can be separated out. 
The question arises: how one could understand the case when a number of remaining modes 
is more then 1 ? 

One point of view is that extra degrees of freedom are inherent to Fokker-type system 
considered literally as a physical model. If complex frequencies are present, the system is 
unstable and thus it cannot form bound states, at least in the AGO approximation. It has 
been pointed out in Appendix D that some modes with real frequencies may contribute 
negatively in the energy what is another kind of instability 
should be adjudged as physically inconsistent. 



In these cases the model 



One can adhere another point of view. We considered two-particle Fokker-type action 
integral where 1-Fokkerians correspond to a free-particle system while 2-Fokkerian describes 
particle interaction. A system of two free particles possesses 6 degrees of freedom. The same 
is true for the interacting system describing within the predictive Lagrangian dynamics. If 
one endows a particle interpretation of the Fokker-type system, the latter should possess 
6 physical degrees of freedom. Thus extra degrees of freedom should be considered as 
a mathematical artifact of the Fokker formalism or specific model, and finally should be 
separated out the physical dynamics of the system. Similar situation arises when considering 



the Lorentz-Dirac equation 40|]. It possesses an extra solution describing exponentially 
accelerating particle even if external forces vanish. This solution is commonly discarded as 
nonphysical. 

The question is how to separate physical degrees of freedom out of nonphysical ones ? 
In our case, how to recognize the only dynamical mode of radial excitations ? A kind of 
selection rule can be suggested if 1) there exist some parameter of nonlocality r such that: 

<l>r{^, xi, ±2) — y 6{^)A{xi, ±2), 

r— >-0 

i.e., in the limit r — > the system turns into a predictive Lagrangian system, and 2) if 
this predictive system admits circular orbit solution. As a consequence, all modes of the 
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dynamical matrix reduce to 5 kinematic and 1 dynamical ones while every extra mode 
disappears. It is possible, for example, if corresponding frequency Iwl — )■ (then this mode 
does not contribute in the Hamiltonian; see ( I2.37p . (ID.lSp and (1D.22P ) or l^l — t-oo (then 
it never can be excited). 

If there is no an explicit parameter of nonlocality, it sometimes can be defined dynamically, 
as a function of the angular velocity Q or of the angular momentum J. This is possible 
because is a parameter with respect to the action I^'^^ in fl4.37p as well as the angular 
momentum J or the quantum number i are parameters as to the classical or quantum 
lamiltonian H^'^^ (12.371) . For example, in the Fokker action formulation of electrodynamics 



30| one can put r ~ f = R^l, where f is a particle speed along a circular orbit of the 
radius R. In the domain f <^ 1 there exists only one mode with the frequency which can be 
identified as Ur{J). By continuity this mode can be recognized and selected in the essentially 
relativistic domain f < 1, while other (extra) modes should be discarded as nonphysical. 
Similarly, one can treat other relativistic systems. 

After kinematic modes are suppressed and nonphysical modes are discarded the subse- 
quent treatment of the two-particle Fokker-type system reduces to the effective single-particle 
case with the effective Hamiltonian 

H,s{J, \Ar\) = + H^^\J, \Ar\) = H^'\J) + COriJ)\Ar\^, (4.65) 

where J = \ J\ = Vj^ and \Ar\ = A*Ar. 

At this point one should refer to a symmetry of the original Fokker-type system which 
is Galilei-invariant in a non-relativistic case or Poincare-invariant in a relativistic case. In 
both cases the effective Hamiltonian is understood as the energy of the system in the center- 
of-mass (CM) reference frame. It is a function of J which is meant as the intrinsic angular 
momentum of the system, and A^ which is the amplitude of interparticle radial oscillations. 
In order to have a complete Hamiltonian description of the system one must introduce 
variables (or operators) characterizing the state of the system as a whole, for example, the 
total momentum P and the canonically conjugated CM position Q. Together with J, Ar 
and the function (14.651) these variable unambiguously determine a canonical realization of 
a symmetry (Galileo or Poincare) group, i.e., the complete Hamiltonian description of the 
system. 
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For example, the total energy of a non-relativistic system is 



H = i(mi+m2)~'P' + H,s{J, \Ar\). (4.66) 

In the relativistic case: 

H = + P\ where M = H,s{J, (4.67) 

i.e., the effective Hamiltonian coincides with the total mass M of the system, while the total 
Hamiltonian (14.671) and other generators of the Poincare group are determined in terms 



of M, J, P and Q via the Bakamjian-Thomas 



quantization of BT model is well elaborated 43|, |4^ . But the spectrum of the mass operator 



BT) or equivalent model 4l|, |42]. The 



can be obtained from (14.651) by means of the substitution (12.381) . 
V. CONCLUSION 

In this paper we proposed a quantization scheme for a two-particle Fokker-type system 
on AGO approximation. For generality it is considered a system which is invariant under 
the Aristotle group, the common subgroup of the Galileo and the Poincare groups. In such 
a way both non-relativistic as well as relativistic systems are included into consideration. 
And only at a very final stage the scheme refers to a genuine symmetry of the system. 

It has been proven that the Aristotle-invariant two-particle system admits a planar circu- 
lar motion of particles with arbitrary (in principle) angular velocity Q provided some rather 
general conditions hold. This is done by the usage of a non-inertial uniformly rotating refer- 
ence frame in which circular orbits are search as static (equilibrium) solutions of equations 
of motion. Radii Ra of particle orbits are stated as certain implicit functions of Q = \ 
Then small perturbations of particle motion around equilibrium point are considered. They 



correspond to almost circular particle orbits 



30l | as referred by the inertial observer. 



The action principle of Fokker type for perturbed motion is derived. It leads to a set of 
linear homogeneous integral-differential equation. General properties of this set has been 
studied. 

It is shown by means of a group-theoretical analysis that a possibly wide variety of 
characteristic modes of this set includes a number of modes which are unessential and must be 
suppressed in order to avoid a double count of degrees of freedom. Also, there is a one mode 
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corresponding to radial inter-particle oscillations. Our attitude is that, namely, this mode 
is physically meaningful. All other characteristic modes (if exist) are unstable or physically 
unacceptable and must be suppressed. The corresponding selection rule is suggested. The 
reduced dynamical system is equivalent to some effective single-body problem. It is put by 



means of the Llosa procedure [25(] into the Hamiltonian formulation which then is expanded 
to a two-body problem in accordance to a symmetry of the origrnal Fokker-type system. 
For example, if the original Fokker-type system is Poincare-invariant, the final Hamiltonian 
description is formulated within the Bakamjian- Thomas (BT) or equivalent model j41, 42| 



on the 12- dimensional phase space P. In other words, as a dynamical system, this BT model 
is a predictive subsystem of the original Fokker-type system. A subsequent quantization is 
straightforward. 

The presented scheme for quantization of Fokker-type models may appear useful in the 
nuclear and hadronic physics. The Fokker-type model of Regge trajectories will be presented 
in a forthcoming paper j46|. 
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APPENDIX A. ROTARY INVARIANCE PROPERTIES OF A SINGLE PARTI- 
CLE LAGRANGIAN 

We use the rotary invariance property fl2.ip of the Lagrangian fl2.2p in the infinitesimal 
form: 

X«L = ..,'(x'=|, + L = 0. (A.1) 

Applying the infinitesimal operators = S\X^ to the Lagrangian fl2.7p we express the 
identities (lA.ip in terms of the variables z, z: 



Then the equation (12.100 . the identities (IA.2p and their differential consequences 



d ~ ~ d ~ ~ 

— L = 0, ^X^L = 0, 
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taken on circular orbit, result in the homogeneous linear set of equations: 

Li = 0, Ei^Vt^Eij^z^Lfi = ZiVt^Lfi^ = 0, 

e^■L^ + ej [zkLij + nh^^z^^Ln^) = e^/Li + e^z^Lij + Zi^TLn^ = 0, 

where z = R. This set permits us to express a one part of coefficients fl2.17p Lj, Lij etc. via 
another part of them. Choosing unit coordinate orts as described before eq. fl2.20p we have: 



Lu 

-%L 



R^2 



LiiL 



12 



-^12 -^22 

^L^ 



^11 -^12 



~ji^2 11^1 



It follows from these formulae the equations fl2.20p - p.22p 



so that 



T - - Lt ~ 



^12 + "^-^2 






APPENDIX B. NONRELATIVISTIC PARTICLE IN A POWER-LAW POTEN- 
TIAL 



Let us consider the following Lagrangian of a non-relativistic particle: 



L = — X — a\x\ = — 7 — aa 
2 ' ' 2 ' 



aiy > 0. 



(B.l) 



Eqs. fl2.12p . fl2.13p and f l2.15p in this case lead to the following circular orbit equation and 
integrals of motion: 

= y-R-\ (B.2) 

m 

J = mR'^n, (B.3) 

Combining (1B.2P with (IB.SP and ( 1B.4P yields the functions which determine the circular 
orbit Hamiltonian: 



.J 



u-2 



z/a 



I/+2 



V 



yum J 



(B.5) 
(B.6) 
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of f B.2 ) 



The radial frequency fl2.25p the present case, = v{v+2)^W ^, simphfies with the use 

Ur = Vi^ + 2n. (B.7) 

It follows from this that the circular motion is unstable at u < —2. 

Gathering ( ]B.6p . f l2.37p . ( ]B.7|) . (IB.SP all together and using the quantization rules f l2.38p 
one obtains the energy spectrum: 



i{i + 1) 



I u 2nr + l 

1 + 



+ 1) 



(B.8) 



It coincides exactly with that formula derived in 
the oscillator approximation. 

In the cases of Coulomb and oscillator potentials we have: 

1/ = -1 : 



45| by solving the Schrodinger equation in 



E = - 
v = 2 
E = 



1 



2nr + 1 



2U + nr + l] 



m 



APPENDIX C. SCALAR ARGUMENTS OF FOKKERIANS 

a = = (Sizi - 8222)^ = (Szi - Z2Y 

= z^ + 2{zi ■ zj[l - cos{n^)] - (n, zi, Z2) sm{n^)} , 

{3a = Xa-X= {-TiZa + V^) " {Za - S^S^Za) 

= {-TiZa + Va) ■ {Za - ZaCOs{m) - n{n ■ Za)[l - COs(fi?9)] 

- {-)''nXZaSm{Q,i!})} , a = 1,2; a = 3 - a, 

6 = Xi-X2= [Si{Zi + Vi)] ■ [82(22 + '"2)] = [S{Zi + Vi)] ■ {Z2 + V2) 

= {zi + Vi) ■ {z2 + V2) cos{Q'd) + (n • Zi){n ■ {z2)[l - cos{Q'd)] 
+ (n, Zi + vi, Z2 + V2) sm{^l'&); 



(C.l) 



(C.2) 
(C.3) 



(C.4) 



here Va = ^XZa is a. vector product of f2 and Za] (n, Zi, Z2) = n ■ {zi XZ2); Sa = exp(taQ) 
(a = 1,2); S = SiSj = exp{m) 
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APPENDIX D. HAMILTONIZATION OF A SYSTEM OF NONLOCAL OSCIL- 
LATORS 



Quadratic terms of the action f l4.37p can be presented in the following simple form: 

^^'^ = ^ E/ [dtdt'p\t)D,i{t - t')p\t'), (D.l) 



kl 



where the matrix kernel D{t — t') = \\Dki{t — t')\\ (here the multindeces are used: k,l = 
{a,i), {b,j) etc.) is invariant with respect to time translations and a time reversal: 
□^(t' — t) = D(t — t'). The time- nonlocal linear equations of motion: 

Y,[dt'DM{t-t')p'{t') = (D.2) 
I ^ 

admit a fundamental set of solutions of the form: p'^(t) = e'^(co')e~"^*. Their substitution 
into the equations flD.2p yields the set of algebraic equations: 

J2Dkiiuj)e\uj)=0, (D.3) 
I 

which constitute the eigenvector-eigenvalue problem for the polarization vectors 6*^(0;) 
and characteristic frequencies u. The latter s are determined from the secular equation 
detD(a;) = in terms of the dynamical matrix: D{u) = Jdt D(t)e''^*. (Here the Fourier- 
image D(ci;) is denoted by the same symbol as the prototype D{t) but with different argument 
which might not lead to confusion). Due to non- locality of the problem ( ID. 21) matrix ele- 



ments Dki{uj) are non-polynomial, in general, functions of u. Consequently, solutions of the 
secular equation form, in general, an infinitely large set of complex and/or real characteristic 
frequencies. Thanks to symmetry properties of the dynamical matrix: 

□^(a;) = D{-uj), Dt(a;) = D{uj*) (D.4) 

this set consists of quadruplets {±ijJa, ±i^a, a = 1, 2, ... } (and/or duplets if Ua G M), and 
a general real solution of the equations (ID. 21) is: 

a 

where = e''{ua) and = e^(a;*). 
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Arbitrary complex variables Aa and Ba parameterize a phase space of the system 
(infinitely-dimensional, in general) which may include both the physical as well as non- 
physical degrees of freedom. Complex frequencies cause necessarily difficulties in a physical 



interpretation of the system which we discuss below (also see 39|). Thus, from the whole 
variety of frequencies we choose only real ones: w* = Ua- The general solution (ID.Sp reduces 
in this case to the following one: 



(t) = J2 {^-e^ e-''^"* + A*Ji e'-'^*} , (D.6) 



where summation spreads over those a for which Im Ua = 0. This is implied in the next 
subsection too. 



D.l. Hamiltonian description: real frequencies 

A current problem at this point is to construct the Hamiltonian description for the vari- 
ational principle flD.l|) on the phase subspace P"* C P*^ of solutions flD.6|) parameterized 



by complex variables Ac,. An appropriate guideline which we adopt for this purpose is the 



Hamiltonian formalism for nonlocal Lagrangians proposed by Llosa and coauthors 
Let us define the time-nonlocal lagrangian: 



24-2 



= k E [dt'p'mUt - t')p\t'), (D.7) 

kl 



in term of which the action (1D.1|) takes a usual form J*^^-' = JdtL{t), and the functional 
derivative: 

^'^(^'^' M) = ^ = i$^p'(t)A.(t-0- (D.8) 



I 



Then, following the Ref. 25|], the Hamiltonian structure on the phase space of the time- 



nonlocal system is defined by the symplectic form, i.e., the closed differential 2-form: 

n = Y[dt [ds Iduxit, s)^^4^^^4^V(w) A 5p\t), (D.9) 
J J op\u) 

where s) = ^(sgnt + sgns) = 0{t)e{s) - e{-t)e{-s), 

Sp'^it) denotes a functional differential of p^{t), and "A" denotes the wedge product. In 
turns, an explicit calculation of the symplectic form determines PER of phase variables (in 
our case - of A^'s). 
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It is evidently that the symplectic form (1D.9P is exact, i.e., Q = SQ, where 



Q = J2jdt jdsxit, s)Ek{-s, t; [p])6p\t) 



(D.IO) 



is a 1-form, so called the Liouville form, defined up to arbitrary exact 1-form (i.e., a total 
differential). A calculation of B (rather than Q) is more simple and convenient for our 
purpose. 

The dynamics of a time-nonlocal system in a phase space is determined by the Hamilto- 



nian 



25| 



H = J2jdtJds x(t, s)Ek{-s, t; [p])p\t) - L{t)\t=o. 



(D.ll) 



Upon integration in eqs. (ID.lOp and (ID. lip the following formula is useful: 



oo s 

dtjdsxit,s)E{-s,t)f{t) = j ds j dt E{t - s,t)f{t). 



-oo 



Let us calculate the Liouville form 9. Using flD.6p and f ID.Sp in (ID.lOp yields 



e 



oo s 

/^^ ldtDu{s)p\t-s)5p\t) 



oo s 

\Y.Y. fdsDkiis) /"rft(A„e^e-''^"(*-^) + A;e|,e''^"(*-^)) X 



-OO 



kl Q/3 



X (^e-''^''*ejM;3 + e''^'3*e^M;) (D.12) 
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which then is convenient to present in the matrix form: 



1^1 'J 



kl a/3 



X dt 



kl otl3 



X dsDklis) 



fc/ a/3 



— oo 



X 



Dlki-^a)-Dlk(-^l3) Dlk{-^a)-Dlk{^p) 



(D.13) 



Due to properties (ID.Sp and (ID.4P all items of the sum flD.lSP vanish except those terms 
which both correspond to a = /3 and include anti-diagonal entries of the square matrix in 
the last line of fID.lSp . Residuary terms possess uncertainty 0/0 which can be eliminated by 
a limit transition: 



kl a 



Dlkij^a) k^A A J Dlkj-K) - Afc(-^a) *k x A* 



Aq, — UJa 



^J2^c.{A*JA^-AJA*J, 



:d.14) 



where 



kl 



du 



(D.15) 



In order to calculate the Hamiltonian we first note that the Lagrangian (ID .711 equals to 
zero by virtue of the equations of motion ( ]D.2p . thus the last term of ( ID.lip vanishes. The 
residuary sum in fID.llI) can be obtained from the Liouville form fID.lOP by means of formal 
substitution 5Aa — )■ —iuJaAa- Thus one gets: 



(D.16) 



If Aq > one can redefine polarization vectors — )■ = Aa^^^e^ in (ID.6P so that the 
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Liouville form and the Hamiltonian simplify: 

= -2 Y^^K^^a - ^o.5Al), (D.17) 

a 

H = Y,^MX- (D-18) 

a 

Then one gets from fID.lTp the symplectic form: 

= 59 = i AM„ (D.19) 

a 

which generates the following PBR: {Aa,Ap] = {A*^,A*^} = 0, {Aa,A*^} = -i6ap. Upon 
quantization Aa — > A^, A^^ — )■ Aj^ one obtains standard annihilation and creation opera- 
tors: [Aq,,^!^] = 6ai3', the Hamiltonian (ID.ISP takes the standard oscillator form and leads 
immediately to the discrete spectrum E = J^a^ai^a + = 0, 1, . . . . We will refer to 

variables Aa, v4* , and the form (ID. 191) as canonical ones. 

Let us now {a} = {a'}U{a"} such that A^' > and A^" < 0. Redefining the polarization 
vectors — )■ = {Aal^^^'^e^ and then the canonical variables Aa" — )■ Ca" = A*^„ reduce 
the Liouville form (ID.14P to the canonical one: 



= -2 J2{A*a,6Aa'-Aa>6A*^,) - i ^ ( A;„M,» - M^,, ) 
a' a" 

= I J2^Al,6Aa'-Aa>5Al,) + i ^(C:„5a»-C„.5C:„). 

a' a" 

The Hamiltonian in this case, 

H = 'Y^Ua'\Aa'\'^ — Ua"\Ca"\'^, 
a' a" 

is not positively defined which feature is characteristic of higher derivative and time-nonlocal 
systems. [39 1. 



D.2. Hamiltonian description: complex frequencies 



If Im cJq, 7^ 0, the corresponding terms in f ID.Sp are unbounded which fact contradicts the 
taken approximation of small p^. One can chose a dumping solution by putting Ba = in 
(ID.Sp . This case however cannot be turned into the Hamiltonian formalism, using the scheme 



by Llosa et.al. 



25| . To see this we consider briefly a general case of complex frequencies. 
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Substituting the solution (1D.5|) into (ID.lOp and accomplishing similar to eqs. ( ]D.12p 



kl 

-^I'^^k x-k . £k _ 



flD.14p (but more cumbersome) calculations yields the Liouville form: 

e = i 5^ {A^iB*JA^-AJB:) + Al{A*JB^-BJAl)} , (D.20) 

a 

with complex (contrary to (ID. 151) ) coefficients 

h 

Then the redefinition ^ = A«'/^e^, ^ = A~^''' in flD3|) reduce (E^O]) to 
the form: 

e = i J] (5:M,-A,55: + A*JB^-BJAl) . (D.21) 
Similarly one obtains the Hamiltonian: 

H = Y, i^cBlA^ + ulAlB^) . (D.22) 

a 

The Liouville form (ID.2ip is not split in variables A^ and B^ which thus are not canonical 
nor appropriate for quantization. Properties of the Hamiltonian (ID.22P in these variables 
are obscured. Thus we change variables into canonical ones. For a brevity we consider only 
one mode corresponding to some quadruplet of characteristic frequencies. Thus hereinafter 
the indices a and summation over a are omitted. 

A choice of canonical variables is not unique. One may choose for one a complex variables 
a, 6 as follows: 

A= {a-h*)/^/2, B = {a + b*)/V2, 

in terms of which the Liouville form indeed becomes separated and canonical: 

e = ^{B*6A - A6B* + A*6B - B6A*) = ^{a*6a - a6a* + b*6b - b6b*). 

On the contrary, the Hamiltonian does not split in a and b modes, 

H = ^{ujB*A + uj*A*B) = ^ {Reuj{a*a - b*b) + ilmuj{ba - b*a*)} . 

In this case, real canonical variables are more appropriate for a quantization. If one chooses: 

A= i-P + ip)/V2, B = {q + iQ)/V2, 

the Liouville form standardizes: B == p6q+P6Q (the notation "md" means "up to a total 
differential" ) while if = —Re a; {qP — Qp) — i Im u{qp + QP) . In this form the Hamiltonian 



has been quantized in [39| and its spectrum is shown to be continuous and unbounded both 



from below and from above. A physical meaning of such a system is doubtful. 
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APPENDIX E. ARISTOTLE-INVARIANCE PROPERTIES OF THE DYNAMI- 
CAL MATRIX 



e,,'fi%,' = 0, (E.l) 



L„,,- + sJ^tL^u - 0, (E.2) 



alj 



L,r, + e,^'n'L^i^^ = 0, (E.3) 

e,j!{zlLai + nh.^z^Un} = 0, (E.4) 

eik'UaLalj + n'siJi'z^^Lanj} + E./lal + eJsij'^n'Lan = 0, (E.5) 

£,,'{4L„,,- + n^^zTLanj} + e,/L^i = 0, (E.6) 

+ Frn^nm + ajn'^-^m} - 0, (E.7) 
Aaij + aJn^A^i- + n,n"{A<,„„,- + aJn'^A-jJ 

+ FrU^-ana.m + sjil'^.t^^m = 0, (E.8) 

+ Pra^.„„i(^l) + e^j,n'%i^^m = 0, (E.9) 

+ Pr? + ^'e^^^a^-am) = 0, (E.IO) 

+ mn'^eJiz'^Aalaj + ^''Si^z'^Aanaj} 

+ 5,/A„,+£,,'£,/l]%^ = 0, (E.ll) 

e^k\<KlJ-rn^e^-z^Aar.J] 

+ nin%gfe'{^^AsUi + JCAanai} 

+ Pr? + ^'e,-z^^-anajm + £,/A„,- = 0, (E.12) 

where Zi — Ri, Z2 — —R2 and Prf = 5f — riin'^. 
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